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Abstract 

A general framework for the solutions of the constraints of pure 
gravity is constructed. It provides with well defined mathematical cri- 
teria to classify their solutions in four classes. Complete families of 
solutions are obtained in some cases. A starting point for the system- 
atic study of the solutions of Einstein gravity is suggested. 



1 Introduction 

This paper is devoted to the search of a general framework for the solutions 
of the Hamiltonian constraints of gravity. The canonical description of Gen- 
eral Relativity is an old problem , which goes on acquiring a renewed and 
increasing interest as the starting point for the establishment of a quantum 
theory of gravity [U [2j [3] , which is probably one of the most ambitious pro- 
grams of modern physics. It is an attempt whose final success seems to be 
still far away. 

The idea to extend the methods of quantum field theories to the case 
of gravity (perturbative approach) came up against a number of difficulties, 
the most evident one being the non renormalizability, because of which the 
number of papers on this topic has diminished in the last years. This is not 
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the case with the non perturbative approach which, on the contrary, maintain 
a considerable appeal. A direct application of the Gauss-Codazzi theorem, 
allows to split the four dimensional spacetime into a 3+1 space and time 
(ADM formalism [U, El [7J , [8], [9]), which provides us with a very useful tool 
to understand General Relativity in terms of our direct experience of space 
and time, i. e. three dimensional slices evolving with a parametric time 
variable. The final result is a canonical theory with a Hamiltonian consisting 
in a linear combination of first class constraints. Then the logical starting 
point is to characterize the manifold where the dynamics is defined. This 
is, precisely, the object of this paper. From the work of these pioneers, this 
problem lives a revival with the introduction of gauge theories of gravity PUJ 
[TTJ [T2J Q3], which are on the basis of the most recent proposals by Ashtekar 
and coworkers. They reach a canonical formalism of gravity, close to standard 
gauge theories, where the fundamental dynamical variables are given by the 
canonically conjugate pair (E?, A ia ). Here E% are the SO (3) densitized triads, 
and A ia the corresponding gauge connections. In the following, latin letters 
a, b, . . . of the beginning of the alphabet are assigned to the coordinates 
defined in the three dimensional slices resulting form a suitable foliation of 
the space-time, while i, j, . . . of the middle of the alphabet are internal SO(3) 
indices (allways in covariant position) running from 1 to 3, as well as a, b, .... 

There have been some attempts to formally solve the constraints [HJ [151 
[TBI ?], thoroughly discussed in a previous reference [T7j, in which the close 
relationship existent between the Ashtekar structure of the constraints and 
the usual ADM one has been emphasized. Using a suitable choice of the 
dynamical degrees of freedom, we get, in this work, a general form of the 
solutions which allows us to classify them in terms of simple and general 
mathematical conditions. As we will see, these imply restrictions on the 
three-dimensional metric. 

With the purpose to offer a complete and unified version of our approach, 
we include, in a first section, a brief review of reference [17], adding besides 
the analysis of some important properties which reveal themselves useful to 
construct a general framework. 
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2 Hamiltonian constraints of General Rela- 
tivity 

The Hamiltonian constraints of General relativity, as derived in the literature 
[20l I2T1 122] , among others, are the weakly vanishing expressions 

V a £« « 0, (1) 
E\F iab « 0, (2) 
e ijk E]E b k F iah « 0, (3) 

describing S0(3) Gauss law, vector and scalar constraints, respectively. For 
simplicity, we have chosen here the Barbero-Immirzi parameter /3 = i, as long 
as at the classical level the Einstein's field equations have the same dynamical 
contents for any value of (3. Ef = ee® are the densitized triads, F iab = 
2d[ a A\i\b] + eijkAj a A kb is the £0(3) field strength tensor, and V a E? = d a E? + 
tijkAjaE^ being the local ST) (3) covariant derivative. Being A ia complex, 
additional reality conditions are required. 

Paying attention to the connections, which are the fundamental dynam- 
ical variables, we have considered in the previous paper |T7] an approach 
which offers an alternative writing of the constraints (DO)-© closer to the 
ordinary geometrical ADM language. To do that, we redefine the SO (3) 
connections as 

Aia Tja ikia j 

where T ia is the part of the connection compatible with the metric and h{ a is 
its intrinsic part, which plays the role of the extrinsic curvature of the three 
dimensional slices. 

With these elements and after some algebra (the reader is referred to 
[T71 \W\ |2"3] [21] for the details), one gets, in coordinate language, the three 



equivalent expressions for the constraints ([I])-© 

k[ab] ~ 0, Gauss law (4) 

D a (k% -5 b Trk) w 0, Vector constraint (5) 

R {3) - Tr(k 2 ) + (Trk) 2 w 0, Scalar constraint , (6) 



where is the scalar curvature defined in the three dimensional space and 
the covariant derivatives are the ordinary Christoffel ones. 

It must be emphasized that ([I])-® and (j3J-([6]) are different versions 
with the same dynamical contents. Three gauges [SO (3), coordinates and 
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reparametrizations] can be arbitrarily chosen to get the final result. The 
problem simplifies drastically after fixing the reparametrization gauge by 
imposing the Dirac's maximal slicing gauge [?] Tr(k) = (other choices are 
possible as we are going to discuss in the last section), so that it reduces to 
solve 

D a k a b w 0, (7) 
R^-Tr(k 2 ) w 0, (8) 

where fc a &is a traceless symmetric matrix. 

It is very easy to verify that the vector constraint ([7]) transforms covari- 
antly in this gauge with respect to the rescaling of g a b and k%, namely 



9ab — > Pgab, 



kt -> v~* /2 k 



what implies 



Furthermore, the vector constraint always admits an identical solution de- 
pending only on the metric tensor. In fact, the Cotton- York tensor C£, 
defined as 

1 



C ab = v acd Dc 



Rj — - S d R 



, „- - (9) 

is symmetric, traceless and identically conserved D a C£ = 0, as required by 
the vector constraint flTJ). 

Consequently, we can establish the general form of k% as the sum of two 
terms 

K = k a h +aCl (10) 

where k% stands for any solution of (j7|) and a is an arbitrary constant. 

This can be easily understood, as far as one looks primarily for solutions 
for k^ in eq. ([7]), which are functionals of an arbitrary metric, they must 
be expressible in terms of the metric tensor and its derivatives up to third 
order. The reason is that the first derivatives of the invariant Tr(fc 2 ), im- 
plicitly present in flTJ), leads through the scalar constraint to a derived scalar 
curvature R^ 3 \ which includes up to second order derivatives of the metric. 

On the other hand, the cancellation of the Cotton- York tensor being the 
necessary and sufficient condition for a manifold to be conformally flat, the 
presence of C£ allows us to distinguish, from the very beginning, conformally 
flat metrics as a very outstanding case, as we will see in the following. 
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3 The choice of the dynamical variables 



A non negligible part of the problems arising in gravity theories is related 
to the search of a suitable set of dynamical variables. In this section, we 
consider a parametrization of k% that reveals itself specially useful to deal 
with the analysis of the possible solutions. For this purpose, we recover the 
triads formalism to write down k% as 



where ky is an ordinary traceless, symmetric matrix. In this way the usual 
5*0(3) symmetry of the triads formalism is recovered. Being symmetric, it 
can be diagonalized with the help of an orthogonal matrix, what enables us 
to deal with its spectral representation. 



where pi are the three eigenvalues verifying p\ + p 2 + P3 = (traceless con- 
dition) and u, v and w are the corresponding eigenvectors. In this way, the 
five degrees of freedom of fc£ arrange as two scalar eigenvalues plus the three 
independent parameters associated with the eigenvectors, which are isomor- 
phic to an SO (3) transformation. The last property will be important in 
what follows. With these assumptions, the final form of k% reads 



where e\ = e°;ui, e?, = e°;vi and e% = e°;wi. 

We impose now reality conditions on the eigenvalues, which become con- 
strained to the real roots of the cubic canonical equation. This restricts the 
discriminant to be zero or negative. A very convenient parametrization for 
both cases is the following 



b ~ e i Kij e jb -i 



hj = UiPiUj + Vip 2 Vj + WiPsWj , 



kl = el pi ei b + e a 2 p 2 e 2b + e^p 3 e 3b , 



(11) 



pi = X cosu;, 

p 2 = — ^-(coscj — a/3 sinu) 



(12) 




where sincu = and sincj 7^ correspond to null and negative discriminant, 
respectively. Now, we have the essential elements to attempt a classification 
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of the possible solutions, paying mainly attention to the vector constraint. It 
is worth to emphasize that as a consequence of such a simple parametrization 
as (j!2p we have formally solved from the very beginning the scalar constraint. 
In fact, a straightforward calculation shows that Tr(/c 2 ), which is independent 
of uj has the value Tr(/c 2 ) = 3A 2 /2. Therefore, taking into account the scalar 
constraint, we have 

^A 2 = i?( 3 ). (13) 



4 Classifying the solutions 

Since there are many and very different solutions of the constraints, a way 
must be found to classify them in a general scheme. For this purpose, we 
rewrite (1TT1) in a more compact notation as 

k% = PlJ e1e jb , (14) 

where pij is the diagonal matrix of the eigenvalues. The vector constraint 
becomes then 

D a (kl) = D a {p ij e a i e jh ) = %ej b d a pij + Pije1D a e jb + pije jb D a e^ = , 
which, after multiplication by e k gives 

e1d a p ik - Pij^/ikj + Pikljij = , (15) 
where the symbols of anholonomy 7^ are given as 

lijk = -e)elD a e ib = e ib e a k D a e) 

and verify 

%k + Ijik = 0, and jkjk = D a e°; . 

Note that the covariant derivative are here the Christoffel ones acting on the 
coordinate indices. 

It is convenient for the later discussion to develop eq. ( 1T51) . which leads 

to 

e a id a pi + 72i 2 (pi - P2) + 7313 (Pi - P3) = , (16) 

e!j0oP2 + 7121 (P2 - Pi) + 7323(P2 - Pa) = , (17) 
eldapZ + 7131 (P3 - Pi) + 7232(P3 - P%) = . (18) 
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Notice that the triads e\ are defined in terms of the matrix elements of k%. 
This implies that the ST) (3) gauge has been fixed in (|IB]) - (118D by imposing 
e", e-2, e-1 to be the eigenvectors of k% with respect to the metric tensor. 

The value of the dynamical variable u in the parametrization ffl2|) pro- 
vides us with a well defined mathematical criterion to classify all the real 
solutions. As will be seen, there are four different classes. 

• Class A. It corresponds to sinu; = 0. This is precisely the case where 
the discriminant of the secular equation cancels. Starting from (Fl6|) - 
(I18p a bit of algebra leads to the following system. 

eI«9 a A + ^A(7 212 + 73i 3 ) = 0, (19) 

e^ a A + 3A 7 m = 0, (20) 
e%d a \ + 3A7i3i = . (21) 

We recall that sin uj = is not the only value leading to the last equa- 
tions, in fact taking sin lu = one recovers (j!9p - (12ip simply by 
redefining the eingenvalues. So that a general definition of this class 
can be formulated by imposing (Tr/c 3 ) 2 = (Tr/c 2 ) 3 /6. 

• Class B. The second class corresponds to cosu; = 0. In this case k% is 
a singular matrix (i. e. det(k) = 0), verifying the equations 

A(72i2-73i 3 ) = 0, (22) 
e a 2 d a \ + A( 7 i2i + 27323) = 0, (23) 
e^ a A + A( 7 i3i + 27232) = 0. (24) 

As in the previous case one remains in the same class when cosu = 
±v3/2, leading to the general definition of this class by means the 
condition Tr k 3 = . 

• Class C. The third class is the general case, in which the characteristic 
equations can be written as 



el«9 a /i+|(3- 5) 7 2i2 + |(3 + 5)7313 = 0, (25) 

e a 2 dMl- 5)]+ M3 - 5)7i2i - 2/^7323 = 0, (26) 

e a 3 d a [v(l + $)] + /i(3 + 5) 7l3 i + 2//5723 2 = 0, (27) 

with 5 = v3 tancu and ji = Acoscj. 
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• Class D. Finally a fourth class occurs when we take directly k% = 
in All H . Therefore, all the solutions in this class depend only on the 
metric tensor. In this case, the distinction between spaces that are 
conformally flat and those which are not acquires a specially relevant 
role, which will be analyzed in the following. 



5 The choice of the coordinates 

Once characterized the different classes of solutions, we can use 3-diff invari- 
ance to choose coordinates. According to a Gauss theorem, this can be done 
by fixing the values of the elements of the metric tensor to render it diago- 
nal. Although other choices are clearly possible, some of the known relevant 
solutions can be written in diagonal form, which on the other hand greatly 
simplifies the calculations. We start, therefore, by writing the metric tensor 
in the form 

/ a \ 

g ab = b . (28) 
\ c / 

We introduce now a "natural" triad associated with this form and constructed 
with the three vectors e\ a = (\/a, 0, 0), e 2a = (0, Vb, 0) and e 3a = (0, 0, y^), 
which are an orthonormal basis with respect to g a \j. Any other basis can 
differ from this one only by the application of an orthogonal matrix Oij , so 
that we can write 

which expresses the relation between the eigenvectors of k% and our "natural" 
triad ej a . Equation (129]) allows thus to parametrize all the arbitrariness of k£ 
simply in terms of a rotation matrix. Inserting (1291) in f|T9|) - fr27|) . it is easy 
to find the final form of the vector constraint in each class. In this way, once 
a metric is adopted, the presence of describes all the generality of the 
theory. 

As a first simple approach, let us consider the case = 5ij. A bit of 
algebra suffices then to find the expression of the vector constraint in the 
four classes. In the case of the class A, it adopts the very simple integrable 
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form 

9! log [\(bcf A ] = 0, (30) 

<9 2 log [Aa 3 / 2 ] =0, (31) 

<9 3 log[Aa 3 / 2 ] = 0. (32) 

As is seen, this implies some restrictions on the fundamental structure of the 
matrix elements of the metric tensor, Aa 3//2 being a function of only the first 
coordinate and A(oc) 3 / 4 = f(x 2 ,x 3 ), f being an arbitrary function of two 
variables. 

The equations in class B can be expressed also in directly integrable form, 
namely 

0, (33) 

0, (34) 
0, (35) 

so that the quotient c/b is independent of the first coordinate while Xc^fa 
and A by/a do not depend on the second and third coordinates, respectively. 
It must be stressed that that, in both cases A and B, these conditions are 
a condensed manner to express families of solutions. Notice, besides, that 
since A = a/2/3 in Dirac's gauge, they involve only the elements of the 
metric tensor. 

The study of the third class C leads us to the more complex system 

<9i log [/i(6c) 3/4 ] -5d 1 log [(o/c) 1/4 j = 0, (36) 
d 2 log [Ml " ^)« 3/2 ] +7^< 9 2 log [a/c] = 0, (37) 

<9 3 log[/i(l + o> 3/2 ] - -A- ^ log [a/6] = 0. (38) 

It is clear that this system is, by no means, so easy to solve as are the previous 
ones. Nevertheless, assuming in general that b ^ c, 5 can be obtained from 
the first equation so that, after substitution in the other two, two conditions 
can be deduced which involve the elements of the metric tensor and their 
first derivatives and define families of solutions. We recall that this results 
correspond to the most simple choice = 5ij , so that this procedure puts 
in evidence the richness of solutions of the problem. 



d\ log \fcjb 

d 2 log [A c^/a] 
d 3 log [A 6a/o] 
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6 The search for the solutions 

In order to find the solutions, we must handle a problem parametrized in 
terms of eight degrees of freedom. Two of them are the eigenvalues of k% 
although, as emphasized before, one is formally found by the condition = 
3A 2 /2. After fixing the gauge, we have three more which correspond to the 
independent elements of the metric tensor. Finally, the remaining three are 
parametrized by the defining elements of an orthogonal matrix, the Eulerian 
angles for instance. We will obtain in such a way the four degrees of freedom 
of pure gravity by solving the scalar and the vector constraints. 

From the mathematical point of view, different possibilities are open. 
Nevertheless, it seems convenient to remain close to our experience of the 
world, so that the best approach is probably to use geometry as a primary 
input. We will consider, therefore, in this section several interesting cases in 
order to test our treatment in geometrical language. 

No doubt, isotropic spaces are obvious candidates for that. Moreover, as 
will be seen later, they allow us to propose an interesting slight modification 
of our approach. An isotropic space is described by a diagonal metric tensor 
that can be written in spherical coordinates as 

0n = l//( r )» 922 = r 2 , #33 = r 2 sin 2 9 , 

f(r) being a function depending only on the radial coordinate r. 
It is easy to see that the scalar curvature adopts the form 

RW = ld r {r[f(r)-l]}. (39) 

Thanks to our parametrization, we can thus write immediately the scalar 
constraint as 

l^ = ^d r {r{f(r)-l]}, (40) 

which expresses a differential relation between / and A. At the same time, 
this defines A as a function depending only on r, a property that highly 
simplifies the calculation, as will be seen. As in last section, we may start by 
choosing the matrix = 5ij, obtaining in this way simple relations enabling 
us to construct solutions of Einstein's equations. For isotropic spaces, a 
simple calculation shows that the vector constraint can be written as the 
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system 



d r (r A cos a;) 



A 



dg COS UJ — 



dg(smuj sin z 



sin 2 fl 

A cLfcosa; + v^3 sinu; 



0, (a) 
0, (b) 
0, fc) 



In the case of sinu = (class A), fHTb) reduces to 
d(r 3 A) = , so that A = 2k Q /r 3 . 
From this and the scalar constraint, it follows 



d r <j r[f(r) - 1] + M 







what leads to the solution 



(41) 



It is straightforward to check that the conditions ( I30l) - (l32j) are clearly satis- 
fied. 

The class D occurs when A = 0. Since the Cotton- York tensor C£ cancels 
for isotropic metrics, only the scalar constraint remains, and it is 



#(3) = -d r {r[f(r) - 1]} = 0, 



what gives the Schwarzshild solution 



f=l- 



2m 



This function describes a static situation since = in this case. 

In the class C, both sinu and cosu are different from zero, so that we 



write (141b ) and (141b ) in the form 



d r uj = cot uj d r log(r 3 A) , dg uj 



-2^3 cot 9 
1 + y/3 cot UJ 
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The integrability condition ded r oj = d r deuo = leads to tancj = — 2y / 3, 
from which some algebra shows that it appears, curiously, an obstruction to 
the integrability as the condition cos# = on the coordinates. The same 
happens in class B. 

This obstruction cos^ = in classes B and C requires a comment. One 
can verify that this is not a property related with isotropy, but rather a 
consequence of the choice of the matrix 0y as a Kronecker 5^, what puts in 
evidence the relevance of the role of the matrix Oij in the construction of the 
different solutions. In fact, there are surely families of solutions corresponding 
to other choices of the matrix. To understand this, one can take Oij as a 
rotation with Eulerian angles (0, 0, ip) such that cos 9 = 0, sin-0 = cos-0 = 
l/-\/2 and arbitrary value of the azimuthal angle, i. e, 

(sin — cos \ 

COS0/V/2 sin0/v/2 l/y/2 . (42) 
— cos 

The result, not detailed here for simplicity, is that the Eulerian angle 
is independent on the azimuthal coordinate tp, what does not pose any 
restriction on the coordinates. 

Isotropic spaces have a property suitable to be used in a slightly different 
approach, due to the fact that the Ricci tensor becomes itself diagonal in these 
spaces. This suggests to parametrize the solution of the vector constraint in 
the alternative form 

k a h = V acd D c M M , 

where the tensor M^b must be symmetric to account for the property of 
kl of being traceless. Moreover, it is easy to show that M db and the Ricci 
tensor R a b commute, otherwise k% would not be divergenceless. Finally, the 
symmetry of k ab requires that 

D a [M2-5 a b Tr(M)} = 0, 

as is easy to see, just cancelling its antisymmetric part. From this, it is 
very easy to obtain again the arbitrariness related to the presence of the 
Cotton- York tensor. In fact, if we take simply 

Ml - 5 5 a Tr(M) = G%, 

Gl being the Einstein's tensor, it turns out that 

M a h =Rl-\ 5 a b R 
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reproduces an identical solution. In this way one can work with a symmetric 
matrix M a & that can be simultaneously diagonalized with the Ricci tensor, a 
property very useful in some cases. 

To end this section, we include a brief example of the inverse problem, 
starting from a well known solution and reconstructing from it the different 
elements of the method. The example is a stationary four dimensional metric 
with axial symmetry and Papapetrou's structure, given as 



9af3 



( 9oo 5f 03 \ 

gn = a 

g 22 = b 

V #03 933= C / 



the components of which are functions depending only on r and 9. The 
extrinsic curvature is in this case 



1 



ab 



V~1oo 



1 ab 



C*i 

a 2 
ol\ a 2 



Solving now the eigenvalues equation with respect to the three-dimensional 
restriction of the metric 

kab^i = Pi9ab^i , 



we obtain the eigenvalues 
Pi = 0, 



Pi = ~P3 = 



P 



'abc v abc 

The corresponding eigenvectors can be easily deduced. They are 




b oti 
2a p 




a «2 1 



2b p ' ^2~ c 



The matrix relating them with the "natural" triad e\ = (1/y/a, 0,0), 
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(0,1/ Vb, 0),e§ = (0,0,l/v^),is 



/ y/aa 2 
P 

b Q(i 

27 



V 



in which one recognizes eq. (l4"2l) if sin i 
is therefore a typical B-class solution. 



P 
a «2 
2_7 

a «2 




1 

i 



2 p v/2 / 
a/o a 2 / p and cos « 



ybax/p. It 



7 Gauge dependence and concluding remarks 

The main result of this paper is the establishment of a general framework for 
the study of the solutions of Einstein's equations in the case of pure gravity, 
which gives a classification based on simple and well defined mathematical 
grounds. For the sake of simplicity, we have presented our results in Dirac's 
"maximal slicing gauge" , in which the problem becomes particularly simple. 
Nevertheless, as is well known, the symmetries of the gravitational Hamilto- 
nian are non internal. This implies that the different gauges maybe more or 
less adapted to the global geometrical properties, essentially open or closed 
spaces. As a consequence, a number of gauge conditions have been proposed 
(ADM, Dirac or York gauges for instance) [2] and [T8] . 

It is not the aim of this paper to discuss the deep dynamical meaning of a 
non internal symmetry, we limit ourselves here to look for a parametrization 
of the problem, in terms of triads language, suitable for the study, from a 
purely mathematical point of view, of the characterization of the possible 
solutions. 

At this purpose, we consider the constraints (4), (5) and (6) as a system 
of equations providing us the natural way to obtain the extrinsic curvature, 
using as input a given triad or metric tensor. This is the reason to adopt the 
basis in which K% is diagonal, in such a way that the problem reduces to the 
calculation of the three extrinsic curvature eigenvalues. 

If we maintain TiK = A different from zero, the parametrization (12) 



15 



and 



reads 

1 A ~ 1 A 

Pi = Acos^ + - A = pi + - A, 

p 2 = ~ (cosco- - v^sintu) + i A = p 2 + ^ A, (43) 

A r . 1 1 

p 3 = (coscj + V3sma;) + - A = p 3 + -A. 

A brief calculation shows that the general form of the constraints (5) and (6) 
becomes 

2 

e^dapi + 72i 2 (pi - P2) + 73i3(Pi - Pa) = g <9 a A, 

2 

e2^ a p2 + 7121 (/Sa - Pi) + 7323(p2 - Pa) = g d a A, (44) 

2 

e39 a p3 + 7131 (p3 - Pi) + 7232 (p3 ~ P 2 ) = - ^3 <9 a A, 

^(3) _ 3 A 2 + 2 A2 = (45) 
2 3 

We notice that, in this parametrization, the scalar constraint is indepen- 
dent on u, thus it can be used to express for instance A as a function of the 
three dimensional scalar curvature and the trace A. In this way, the problem 
reduces to discuss the vector constraint (jHJ) as defined by a set of partial 
differential equations depending on two unknown functions, A and to. There- 
fore, to solve them an additional condition (gauge) in A, uj or both is needed. 
We emphasize that, the eigenvalues being scalar quantities, a condition of this 
kind is invariant with respect to general coordinate transformations in the 
three-dimensional space. 

The structure of the vector constraint strongly suggests to fix the gauge 
giving a condition on the trace A. In this manner, the vector constraint 
becomes an ordinary differential equation for u. As a matter of fact, the most 
common choices occurs precisely when A is taken as a constant function, here 
of course "constant" means a function independent on the three dimensional 
space coordinates, a true constant for instance or a function depending on 
the time variable (Dirac or York proposals [22] )• m both cases the second 
term in eq. (1441) vanishes, so that we recover eqs. (I16j) - (jl8j) and therefore 
our classification holds. It must be emphasized that, the scalar constraint 
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being independent on uj, it is "formally" solved, once the value of A is fixed 
as long as A is algebraically given in terms of A and the three dimensional 
scalar curvatures. 

Therefore, it remains only the vector constraint as the natural criterion 
to classify the different solutions of the problem, a classification that, in our 
scheme, is simply given in terms of the possible roots of the cubic secular 
equation. A criterion on the other hand that remains valid in any possible 
gauge. In fact, it must be emphsized that eq. (j4"4"|) can alternatively be 
used to directly characterize the different classes. Not only that but, due to 
the fact that A and B cases correspond to fixed values of u, the problem is 
formally solved without any additional condition. However, Dirac's gauge is 
an useful frame to illustrate the richness of solutions as long as it describes, 
in a natural way, the relevant role of the Cotton- York tensor in the search 
for three dimensional geometrical configurations. 
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